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^ , Abstract. We study the dynamical stability of the macroscopic quantum oscillations 

^ I characterizing a system of three coupled Bose-Einstein condensates arranged into an 

open-chain geometry. The boson interaction, the hopping amplitude and the central- 
well relative depth are regarded as adjustable parameters. After deriving the stability 
' diagrams of the system, we identify three mechanisms to realize the transition from 

an unstable to stable behavior and analyze specific configurations that, by suitably 
I tuning the model parameters, give rise to macroscopic effects which are expected to 

' be accessible to experimental observation. Also, we pinpoint a system regime that 

O . realizes a Josephson-junction-like effect. In this regime the system configuration do not 

depend on the model interaction parameters, and the population oscillation amplitude 
is related to the condensate-phase difference. This fact makes possible estimating 
the latter quantity, since the measure of the oscillating amplitudes is experimentally 



4 : accessible. 
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PACS numbers: 03.75.Kk,03.65.Sq,74.50.-f r,05.45.-a 



Since the first observation of Bose-Einstein condensates (BECs) in dilute weakly 
interacting gases of bosons in 1995, great efforts both experimental and theoretical 
have been addressed on this subject. Thus nowadays, thanks to the impressive progress 
of experimental techniques, BECs represent a source of inspiration for new challenging 
problems in quantum physics. A notable aspect inherent in the dynamics of coupled 
BECs is that the nonlinear character of their equations of motion entails a particularly 
rich phenomenology where nonlinear effects including chaos can be studied in a quantum 
environment. Among the latter, in recent experiments have been observed phenomena 
such as super- fluidity |T] , Josephson tunneling [2] and atom optics [3] . 

If the recent past has been prolific of experiments on BECs arranged into arrays 
with huge number (of the order of several hundreds) of condensates, at the present, 
the experiments on chains with a few interacting BECs (of the order of several unities) 
has received a limited attention [l]-[8] despite the increasing amount of theoretical work 



Control of unstable macroscopic 



2 



devoted to such systems. Over the last few years, work on small-chain models has been 
focused in particular on the two- well system (dimer) [9]- [13] and the three- well system 
(trimer) [H], [15] while their many- well generalization has been studied in |16j-|21j. 
The theoretical prediction of macroscopic effects such as coherent atomic tunneling [9], 
population self-trapping [22], and Josephson junction effect [23], have been confirmed 
by recent observations in the experimental realization of the dimer [l]-[8]. 

These results stimulate the theoretic study of the BEG small chains in order to 
predict new macroscopic effects that, in addition to furnishing a deeper insight of 
stability properties and an increased control of systems, open the possibility of designing 
significant experiments. In this respect, the recent experimental technique discussed in 
[51 H] and successfully applied to realize a two- well system provides an extremely effective 
tool for engineering small arrays. The superposition of a parabolic confining potential 
to a linear optical lattice allows one to select two wells (dimer) with possibly different 
depths by changing the parabola position and amplitude. Likewise, a simple, suitable 
change of the latter should enable one to select M-well arrays, with M arbitrary, whose 
(local) potentials feature a parabolic profile. The trimer case is thus very promising in 
that, in addition to represent a practicable experimental objective, it exhibits behaviors 
typical of more complex arrays. In fact, while the simplest BEG array (the dimer) 
features an integrable dynamics [21] , in the trimer case the apparently harmless addition 
of a further coupled condensate is sufficient to make the system nonintegrable. As a 
consequence, the trimer displays strong dynamical instabilities in extended regions of 
the phase space [25] and a whole new class of behaviors including chaos. In general, 
unstable behaviors have been found in condensates subject to kicked pulses [36] or to 
an anisotropic trapping potential [37] and within soliton dynamics [38] . 

Goncerning the trimer, a conspicuous number of recent papers has addressed various 
aspects of its dynamics. These include, within a mean-field approach, the occurrence of 
a transistor-like behavior [26], and the self-trapping [27] in a triple- well asymmetric trap 
of repulsive bosons, the ring-trimer dynamics with attractive interaction in the presence 
of a single on-site defect [28] , and the chaos control through an external laser pulse [29] . 
Recently, the dynamical instability has been studied in terms of quasimomentum modes 
for a 5- well semiclassical Bose-Hubbard model including a linear potential |30j. Also, 
more related to purely quantum properties, second- quantized models of the trimer have 
been considered to study the properties of quantum states relevant to interwell particle 
exchange [31] and entanglement [32], the semiclassical quantization of the Bogoliubov 
spectrum [33], the structural changes of trimer eigenstates in terms of the tunneling 
parameter [34], and the occurrence of vortex-like states [35] . 

In the present paper we study the classical dynamics of an asymmetric open trimer 
(AOT) made of three coupled BEGs arranged into an open chain geometry, where both 
the interwell tunneling amplitude T and the relative depth w of the central well, are 
regarded as adjustable parameters. Our objective is to evidence several macroscopic 
effects that can take place in the AOT dynamics to provide a useful guide for future 
experiments. In section |2] we introduce the semiclassical equations of the trimer mean- 
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field dynamics. In section 12.11 after determining tlie solutions relevant to the fixed 
points of motion equations, i. e. the proper modes of trimer dynamics, we focus mostly 
on establishing, via standard procedures, their stability character for experimentally 
significant values of effective parameters r := T/UN and u := w/UN, where U and 
are the effective interatomic scattering and the total boson number, respectively. In 
section [3] we review the (three) stability diagrams (one for each class of fixed points) 
that summarize characters and properties of fixed points. These supply an exhaustive, 
operationally valuable, account of the system stability properties and allow to determine 
the stable/unstable charater of the system evolution for initial conditions situated, in 
the phase space, close to a given fixed point. Both the derivation and the use of such 
diagrams, rapidly sketched in 09], are discussed here (and in the relevant appendices) 
with some detail together with various new applications. 

By using the stability diagrams as a map for detecting the critical behaviors of 
trimer, in section H] we find parameter-tuned macroscopic effects (both in the dimeric 
and nondimeric regime) and evidence some features that may prove experimentally 
relevant. Such effects are shown to arise when crossing either the boundaries of regions 
with different stability characters or the boundaries of forbidden regions or, finally, by 
exploiting the coalescence/bifurcation mechanism characterizing certain fixed points. 
In section [5] we consider the special case of vr-like states relevant to the fixed points 
of the regime charaterized by a central depleted well. We investigate the unexpected 
relation between the average population oscillations and the initial phase difference of 
the lateral condensates, suggesting a possible experiment. Finally, section [6] is devoted 
to make some concluding remarks. 

2. Dynamical equations of the AOT and fixed points 

The quantum Hamiltonian of the AOT on a linear chain reads [T6] 



where site boson operators a], satisfy standard commutators [oj, a{,] = 6ik, Ui = afat, 
V is the energy offset of the wells, w is the central-well relative depth, T is the hopping 
amplitude and U ^ 2TTah'^/m {a the scattering length, m the boson mass) is the on-site 
boson interaction. The total boson number N = rij, commutes with Hamiltonian ([1]). 
The exact form of parameters U and T is calculated in [ID] showing how such quantities 
depend on m, on the laser wavelength A, and on the potential depth Vq of potential 
Vop{x) = Vo sin^(27rx/A) determining the optical lattice. The 3-well linear system, 
obtained by superposing the optical potential to parabolic potential Vp{x) = 
confining bosons [1], [5], is expected to have lateral wells with depth v smaller then the 
central- well depth v + w due to the presence of the parabolic potential. The dynamics of 
the AOT with strong tunnelling amplitudes {T/U » 1) and large average numbers of 
atoms per well can be described by the mean- field form of ([1]). The method developed in 
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combines the coherent-state picture of quantum systems with the time-dependent 
variational approach [32] providing the semiclassical version of quantum model in 
which the dynamical canonical variables are the expectation values Zi := {Z\ai\Z), 
z* := {Z\af\Z). The so-called trial macroscopic wave function \ Z) is standardly assumed 
to be \Z) = Yli where 1 < ^ < 3 and \zj) are the Glauber coherent states such that 
ai\zi) = Zi \zi) for each i. Then, the resulting mean-field Hamiltonian reads 

3 

n, = Y, [^k.f - ^k.f ] - Mz2? - |[^2*(^i + ^3) + c.c] , (2) 

which, through the Poisson brackets {z*,Zk} = ji^jk, yields the dynamical equations 



ihzj = {2U\zj\'^ - v) Zj - ^Z2 j = 1, 3 

ihZ2 = i2U\z2\^ - V-W) Z2-J {Z1 + Z3) 



The validity of such a semiclassical picture is confirmed experimentally [4J, [5J for 
the time scale (hundreds of milliseconds) of the macroscopic oscillations in a two-well 
system. The effect of approximating the system state \Z) through a product of localized 
(coherent) states Yli appears on a much longer time scale [33] when \Z) must be 
expressed as a superposition of Fock states \n) = |ni, n2, na) {afai\n) = nj|n)) to include 
(quantum) delocalization effects. An early analysis of equations ([3]) for w = is reported 
in [33], [35]. It is easy matter to check that system ([3]) embodies the conservation of 
boson total number = namely {N,H} = 0. Constant can be exploited to 

introduce new rescaled variables Zj = Zje~^'"^^^/ Vn whose evolution is governed by 

zZj = 2|Z/Z,-iZ2 J = 1,3 

^ tZ!, = {2\Z2\^-u)Z2-^iZ, + Z,) ^ ^ 

where 

iy = w/{UN), T = T/{UN), i=UNt/h, (5) 

and the prime denotes derivation with respect to the rescaled time variable t. This 
shows that the two significant parameters of the trimer dynamics are u and r (namely, 
at fixed A^, w/U and T/U), parameter v being incorporated in the time dependent 
factor e~™*/^ of new variables Zj. Furthermore, it is easy to verify that the condition 
zi(t) = zs{t) is valid for t > (or, equivalently, that Zi(t) = Z^it) holds for every t > 0) 
if, initially, Zi{0) = z^^O). Such a condition identifies an integrable subregime where the 
lateral condensates have the same population and phase (see, e. g., [33], [36] and [3Z|) 
which in the following we shall refer to as dimeric regime. 



2.1. Fixed points 

We can get insight about the behaviour of the system in different dynamical regimes 
through the study of special exact solutions along with their stability character. With 
the latter we mean the simplest (single mode) periodic orbits in the phase space of the 
form Zj[t') = XjC*™*'. Substituting the latter in equations (j3]) provides the equations for 



Control of unstable macroscopic 



5 



variables Xj. Such equations can be equivalently obtained from iZ'j = = {Zj, H—mN} 
with H = T-C-^{Zi, Zl), where the chemical potential m, a Lagrange multiplier, allows 
one to include the conserved quantity = = 1 derived from the total boson 

number N = {zil"^. Owing to the condition Z'i = we name fixed-point equations 
the equations determining Xi, i = 1,2, 3. These have the form 

T T 

= (2 - m) Xj - -X2, = [2xl-m- u) ^2 - ^ (a^i + 0:3) (6) 

where j = 1,3 and the fixed-point coordinates Xj are normalized to unity: Ylij^"] — ^■ 
In equation [6] we have exploited the system symmetry, inherited from equations HI under 
a global phase shift Xj ^ Xjexp(2$) to restrict coordinates Xj (which are in principle 
complex numbers) to the real axis. Let us summarize the solutions of system ([6]). 

Central depleted well. This parameter-independent solution is obtained when the central 
well is strictly depleted, that is X2 = 0, whereas the lateral ones have the same population 
and opposite phases. Due to the constraint on the number it must necessarily be 
Xi = —Xs = -^^fmpl = ±l/v^, m = 1. In summary 

x, = -xs = ±l/V2, X2 = 0. (7) 

We shall refer to this solution as central depleted well (CDW). Since configurations 
where either one or both of the lateral condensates have zero population are forbidden 
by equations (E]), we are left with the following two possibilities. 

Dimeric solutions. Similar to ([3]) and (jl]), equations ([6]) possess dimeric solutions 
characterized by X2 7^ and xi = X3 7^ 0. These solutions can be expressed via 
the polar representation 

Xi = x^ = cosO/ V2 , a;2 = sin^. (8) 

Introducing this parametrization in ([6]) one finds that the tangent tan^ =: a of the 
angular coordinate 6 identifying each fixed point is given by the real roots of 

V{a; T, u) = a' + —a' - - 1 = , (9) 

T T 

with A,y = v^(2 — v). By = \/2{l + v). In the following we shall refer to this class of 
solutions as dimeric fixed points (DFP) whose derivation is discussed in Appendix A.l 



Non-dimeric solutions. The remaining single- mode solutions of system ([6]), characterized 
by xi 7^ Xs, Xj 7^ 0, will be referred to as non-dimeric fixed points (NFP). In this case 
it is appropriate to switch to the alternative set of coordinates Xi = (xi + X3)/^/2, 
X2 = X2, and X3 = (xi — X3)/v^ before setting 

Xi = Rcose, X2 = Rsm9. (10) 

In this case the radial coordinate R reads 



Bl \ , ;i + a')(2V2-ar) X, 

272(2 + ' " = '-»-x7 (") 
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and depends on the Hamiltonian parameter r and on parameter a. By introducing this 
parametrization in system ([6]) one finds that a is given by the real roots of 

P{a;T,iy) = a^ + —a^-—a-^ = 0, (12) 
r TO 

where = 2^2(2 - z/)/3, = 4^2(1 + z/)/3. Since the physical solutions of system 
([6]) have to be real it is clear that the roots which make R{a) or 



/ X / Z ^ /4+V2«(l + ra2) 

imaginary, must be discarded. The next section contains a qualitative discussion of the 
three stability diagrams describing the entire set of solutions of system ([6]) for any choice 



of T and u. A detailed derivation of such fixed points is given in Appendix A. 2 



3. Stability diagrams 

The stability diagrams introduced in [39] allow one to know both number and stability 
character of the fixed points of trimer dynamics for any choice of significant parameters 
r and u. In the previous section we have shown that vector x = (xi,X2,X3) describing 
a fixed point can be parametrized by the angular coordinate 6 of the two-dimensional 
polar representations ([H]), and ( fTOj) and ( ITTl) for the DFP and the NFP, respectively. To 
determine parameter 6 (and thus the vector x relevant to a given fixed point) we must 
find the roots of the fourth-degree polynomial 

V{a;T,iy) = a^ + —a^ + —a + C = (14) 

T T 

in the unknown quantity a = tan^. This allows one to determine a = X2/V^Xi in the 
dimeric case with C = —1 (see Qj), and a = X2/V^{xi + xs) in the non-dimeric case 
with C = —4/3 (see (fT2|) ). In the latter case vector x is obtained by using the constraint 
Xi+X2 + xl = 1. Therefore, for a given value of parameters r and u, the number of 
roots of f|T^ establishes how many DFP and NFP there exist. This information is 
summarized in figures [H and [3l 

The stability character of the fixed points solving equations is depicted in the 
stability diagrams of figures [2] and H] for the DFP and the NFP case, respectively. In 
such figures any point in the plane 6' — r is implicitly associated to a definite value of 
parameter i/, owing to equation f|T4l) . This can be seen by rewriting f|T^ as 

[A,t&n^{e) + B,] tan(g) 
"t^""^^)' = tan^(^) + C ' ^''^ 

which, for a given value u = i), defines a curve r[tan(0); z>] in the plane 6-t consisting of 
two or more non-connected branches, both in the dimeric (C = —1) and the non-dimeric 
case (C = —4/3). At the operational level, in both figures, for a given choice r = f and 
u = 9, configuration parameter 6 can be obtained by intersecting the curves r[tan(^^); z/], 
for the dimeric and non-dimeric case, and the straight line r = f. Once (triplet r, u, 
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9 of) a fixed point has been identified, its stable/unstable character can be evinced 
from figures [2] or H] thanks to the colour corresponding to a given value of 9 and r (see 
the relevant figure captions). The linear stability analysis determining the character of 
fixed points (and, in practice, the whole structure of such phase diagrams) is discussed 
As an example, in figures [H and [3] are reported with different dashing 
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Figure 1. Number of dimeric fixed points on varying parameters r and v. The border 
of the lobe where four solutions are found is described by (jA.Sp . 
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Figure 2. Stability diagram for the dimeric fixed points (DFP). The fixed points 
belonging to the white, medium grey and dark grey regions are always minima, unstable 
saddles and maxima, respectively. The latter two regions are divided by the curve 
T = d{ta,ii9) defined in [Appendix A.l| The tip of the "parabolic" branch moves 
along the white dashed curve. The white dotted and dashed curves are (ii(tan0) 
and (i2(tan(?), respectively (see [Appendix A.Tj for details). Note that the curves drawn 
with diffrent dashing styles, correspond to the five values of v selecting the straight 
horizontal lines reported with the same dashing styles in figure [1] 

styles, five straight horizontal lines corresponding to several values of the v. The five 
curves drawn in figures [2] and HI refer to the straight lines drawn with the same dashing 
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Figure 4. Stability diagram for the non-dimeric fixed points (NFP). The fixed 
points belonging to the light, medium and dark grey regions are always stable saddles, 
unstable saddles and maxima, respectively. Note that the white solid curves are the 
loci of the tips of the two parabolic branches. The patterned regions must be ignored, 
since the solutions found there yield complex Xj's (see [Appendix for details). The 
curves drawn with diffrent dashing styles, correspond to the five values of v selecting 
the straight horizontal lines reported with the same dashing styles in figure [3l 



styles in figures [T] and [3l for tlie dimeric and non-dimeric case, respectively. Notice that 
in general these curves consist of two or more non-connected branches. In the dimeric 
case there are always two unbounded branches featuring an asymptote at 6* = +7r/4 or 
at 9 = — ± TT / 4. For u G [—1, 2] a further, roughly parabolic, bounded branch appears. 
This means that for any choice of the parameters r and u there are always at least two 
DFPs, and that, for suitable choices of these parameters, two more fixed points appear. 

In the non-dimeric case, figure HI the curves r[tan(6'); z/] consist of either one or 
two bounded parabola-like branches, hence featuring zero, two or four intersection with 
a constant-r straight line. However, for certain values of the parameters, not all of 
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these intersections represent fixed points, i.e. real solutions of system Indeed, as 



it is discussed in Appendix A. 2 some of the real roots of polynomial V{a] r, u) -more 
precisely the ones lying in the patterned regions of figure HI- result in complex- valued 
solutions of system (Q and therefore must be discarded. In particular, the patterned 
(forbidden) regions of figure H] can be shown to correspond to the white region in the t-u 
plane of figure [3] where no solution is permitted. Figures [2] and H] also supply important 
information concerning fixed points and their characteristic energies. Referring to the 
different colours used in such figures, a fixed point is a (stable) minimum, a stable 
saddle, an unstable saddle or a (stable) maximum of the energy function depending on 
whether the region it lies in is white, light grey, medium grey or dark grey. 

As we mentioned above, other than the DFPs and the NFPs, a further fixed 
point, the CDW configuration, is present for any value of r and u. The parameter 
independence of the relevant configuration x = (l/-\/2, 0, — 1/-\/2), makes the use of a 
polar representation unnecessary, and allows one to draw its stability diagram in the 
parameter plane t-v in a. direct way. This is illustrated in figure [5l 




Figure 5. Stability diagram for the CDW saddle point. Note that large stability 
regions appear (light grey filling). The boundaries are the straight line v = Q and the 
curves r — ci{v) and r — C2{v), described by (jA.lip . The asymptote = 1 of the 
latter curve is also displayed (dashed line). 



4. Macroscopic effects in the AOT dynamics 

A careful look at stability diagrams [2] and HI in addition to provide, at first glance, 
number and properties of the fixed points of AOT dynamics for any choice of r and 
z/, allows one to single out at least three ways to produce macroscopic, and hence 
experimentally observable, effects by means of a simple change of model parameters. 

Crossing boundaries of regions with a different character. For suitable choices of 
z/, curves r[tan(5);z/] may cross the boundary dividing regions with different stability 
properties. Thus, by slightly and adiabatically changing the Hamiltonian parameters 
(e.g. the optical potential amplitude) the system can modify significantly its behavior. 
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This is the case of (one of) the asymptotic branches of the dimeric curves, which, for 
any u > 2\/lA — 7 ~ 0.48, crosses the boundary between the maxima and the unstable 
saddles region (see, e. g., the branch in figure [1] corresponding to v = 0.6 in figure [3]). 
Likewise, in the non-dimeric case, one of the parabola-like branches of curves r[tan(^); u] 
always crosses the interesting unstable "isthmus" splitting the large stable (light-grey) 
region in figure HI provided that z/ < 2. 

Coalescence/ bifurcation of fixed points. At the "tips" of the parabola-like branches 
of curves r[tan(6'); z/], a coalescence/ bifurcation effect oi pairs of fixed points takes place. 
Hence, for some values f of the hopping parameter it is possible to find a value i) with 
two coincident solutions of fixed-point equations located at the parabolic-branch tip. By 
varying z/, the branch tip moves away from f, and the coincident solutions either split 
(i.e. bifurcate) or disappear accordingly. Let us emphasize that in the non-dimeric case 
(figure HI) the coalescence phenomena always involve fixed points with different stability 
characters, whereas, in the dimeric case (figure [2]) only the coalescence of two fixed 
points of the same kind may happen, provided that u > 2\/l4 — 7 (see next section). 

Crossing forbidden-region boundaries. Significant changes in the phase-space 
structure of the system are expected as well where the non-dimeric curves r[tan(^); u] 
cross the boundaries of forbidden zones. Indeed, in this slight change in the 

parameters may cause the sudden appearance/disappearance of a fixed point. As to the 
CDW fixed points, which -according to the Hessian signature- is always a saddle, figureO 
clearly shows that it may display an either stable or unstable character, depending on 
the choice of parameters r and u 

Based on these three mechanisms, we discuss the results of several numeric 
simulations evidencing how macroscopic effects are induced in the AOT dynamics by 
slight changes in parameters r and u (similar effects on the stability of the many-well 
system with U < 0, of vortexlike modes in a ID array and of localized modes in a 2D 
array have been studied in [18], [20] and [21], respectively). These entail that trajectories 
based on the same initial configuration may exhibit an either regular or chaotic behavior. 
The character of the AOT dynamics, can be conveniently investigated by constructing 
appropriate Poincare sections within the four- dimensional reduced phase-space spanned 
by the (reduced) set of canonical variables 0i,^i,02,^2- This reduction is obtained via 
the canonical transformation {zi,z*) {(pi, (j)2,4',^i,^2, N), where 

|2)/Ar, 



01 = ^2-^1, 


6 = 


(k2| 


2 + 


l^sl — 


\zi\' 


02 = ^3 — ^2, 


6 = 




2 _ 






^ = es + 9i, 


= 




l' + 


\Z2? + 


^31 



)/N, (16) 



once the condition A^ = const is taken into account and one notices that the ip evolution 
is nonautonomous (see for details on the reduced dynamics). 

We have numerically integrated the trimer equations of motion in the reduced 
coordinates, by a first-order bilateral symplectic algorithm. Moreover, we have 
qualitatively investigated the degree of regularity of the ensuing trajectories by means of 
Poincare sections. The latter are obtained by projecting the points belonging to a given 
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submanifold, characterized by a given value of one of the four reduced variables, onto 
the plane spanned by two of the three remaining variables. All of the Poincare sections 
reported in the following refer to a set of trajectories based on an array of starting points 
surrounding a significant configuration of the system, usually a fixed point, of reduced 
coordinates {4>i,^i,(f>2,^2)- In every instance the points belonging to the submanifold 
^2 = have been projected onto the (f>i-^i plane. 




0.4 0.8^1.2 1.6 2.0 



Figure 6. Parameter-induced macroscopic change in the stabihty of the trimer 
dynamics. The figure shows a portion of the lower half plane of dimeric stability 
diagram [21 along with the dimeric r-curves (see (HSJ, (|A.6|) ) relevant to = j^m = 2/3 
(solid) and i/ = 1/3 = 17/24 (dashed). The white dotted and dashed curves appearing 
in figure [2] have been omitted for clarity. From the figure it is clear that, for both 
of the considered values of v, there exists a suitable value of the hopping parameter 
— tm = 7/4-\/6 and ts = 13/8\/6, respectively — producing a solution a,t 9 = —tt/3 
(horizontal dashed line). Hence slight differences in the parameters may result in 
different stability characters for the same configuration {9 ~ — 7r/3). Indeed M is a 
(stable) maximum, whereas S is an unstable saddle. 

4.1. Priming instability in breather-like regular oscillations 

The first macroscopic effect we discuss here has been briefly discussed in |39] and 
concerns the dimeric configuration x = (l/-\/8, — l/v^S)*, -with bosons in 
each lateral well and the remaining 3N/A in the central one- which is represented by 
the angular variable, 6 = —ir/S. One can show that the latter is a (dimeric) fixed 
point provided that 4:^/6u + 8r — 5\/6 = 0. Furthermore, such a fixed point is either a 
maximum or an unstable saddle depending on whether u is smaller or larger than the 
critical value 9\/6/32. This is clearly visible in figure [6l where two different parameter 
choices (tm, i^m) and (rs, z/g), shown with the letters M and S, are considered. 

The crossing of the boundary separating the stable and unstable regions has 
evident consequences. The Poincare sections (reported in |39]) of trajectories, based 
on a set of points closely surrounding the fixed point under concern (^i, ^1, 02, C2) = 
(tt, 3/4, — vr, — 3/4), mirror the dramatic change in the dynamical behaviour involved 
in such a slight parameter change. While choice {tm,i^) results in the neat and 
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regular Poincare sections characterizing a stable fixed point, with choice (rs, z/s) sections 
densely fill an area considerably larger than the previous example, despite the relevant 
trajectories are based on the very same points [39] . 

4-2. Loss of anti-breather coherence in the isthmus domain 

An instability outbreak similar to that discussed above is found in the neighborhood of 
the "instability isthmus" of the non-dimeric stability diagram |H Indeed, the possibly 
present upper branch of the non-dimeric r curve always crosses such unstable region. 
This is the case of four of the five curves plotted in figure H] and of the curve with = 
plotted in figure [3 The inset of figure [7| is a magnification showing fixed points A, B and 
C, located in the proximity of the instability "isthmus" for three slightly different values 
of r (ta = 0.22, tb = 0.29 and rc = 0.36). These fixed points correspond to (slightly) 
different values of coordinate 9, thus the corresponding configurations are quite similar 
to each other. In these three cases more than 97% of the total population is almost 
equally shared by the central and one of the lateral condensates. The phases of these 
condensates are parallel to each other and opposite to that of the remaining, almost 
empty, condensate. This kind of configuration corresponds to the neighborhood of the 
reduced-phase-space point with ^i, 02, ^2) = (0, 0.016, vr, —0.984). Figure [8] shows 
the Poincare sections based on a set of points closely surrounding it, for the three choices 
of r listed above, ta, tb and rc- The first and last choice -corresponding to trajectories 

TT 
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_ TI 

4 

_ TI 

2 1 r 3 4 

Figure 7. Parameter-induced macroscopic change in the stability of the trimer 
dynamics. The figure shows the magnification of a portion of the upper half plane of the 
non-dimeric stability diagram, figure [4l along with the non-dimeric r-curve (see (jl5l) . 
(|A.9P ) relevant to = (solid) and the straight lines corresponding to three different 
choices of the hopping parameter, ta, tb and tq (dashed). The ensuing intersections 
are denoted A, B and C, respectively. Notice that B belongs to the narrow unstable 
region (which we term isthmus)^ whereas A and C belong to the wide stable region 
surrounding it. The changes in the stability character of the fixed point resulting from 
such a relatively slight increase of r are evident in figure [H 

based in the neighborhood of a stable fixed point- are both characterized by extremely 
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Figure 8. Poincare sections of trajectories based in the vicinity of the configurations 
A, B, C selected in figure [T) In every case one of the lateral condensates is almost 
empty, whereas the other two share almost equally the remaining boson population, 
which adds up to more than 97% of the total. The phases of the macroscopically 
filled condensates are parallel to each other and opposite to that of the almost 
empty condensate [such configurations are close to the reduced-phase-space point 
^1, 02, C2) = (0, 0.0f6, TT, —0.984)]. The three situations are distinguished by (from 
left to right) r = 0.22, r — 0.29, r — 0.36; in every instance ly = 0. The features of 
the Poincare sections clearly mirror the stability character of the relevant fixed points. 
Notice furthermore that the central Poincare section, relevant to the choice falling in 
the unstable "isthmus" of figure [71 covers a region roughly an order of magnitude wider 
than the lateral sections, relevant to stable regions. 

regular sections where the populations "flicker" about their initial configuration. In the 
remaining case, the fixed point belongs to the "unstable isthmus" (see figure El), and the 
consequent irregular trajectories produce a quite different pattern. Also in the present 
case, a slight change in the parameters destroys the steady, flickering populations. Once 
again the Poincare sections are indistinguishable from one another and densely flll an 
area considerably larger than in the previous cases. 

4.3. Instability suppression due to fixed-point disappearance 

So far, we have examined two situations where the change in the phase-space structure 
is associated with a modiflcation in the stability character of a flxed point. Nevertheless, 
such a change may be also related to the disappearance of a fixed point. When z/ = 0.52, 
e. g., a consistent portion of the lower branch of the non-dimeric r-curve falls in the 
forbidden zone, similar to the curve featuring the narrowest dashing style in figure HI 
However, for r = 1.66 (i.e. very close to the branch tip), both of the solutions associated 
to the lower branch are acceptable. According to the stability diagram, one of them is 
a maximum, whereas the other is an unstable saddle. Actually, the trajectories in the 
neighborhood of this configuration, produce the expected irregular Poincare sections, 
displayed in the central panel of figure M As r is decreased, the unstable solution moves 
toward the forbidden zone (see figure Hj). For r = 1.5, the maximum is the only non- 
dimeric fixed point surviving. The change in the phase-space structure associated to the 
disappearance of the unstable fixed point is quite evident in the leftmost panel of figure El 
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Figure 9. Regularization of the phase-space due to the disappearance of an unstable 
fixed point. All of the three sets of Poincare sections refer to = 0.52, and have been 
obtained from trajectories based in the vicinity of the same configuration. For r = 1.66 
such configuration is an unstable saddle, and produces quite irregular trajectories, as 
shown in the central panel. Slight variations of the hopping amplitude cause the 
disappearance of the unstable fixed point, and this in turn results in an increase in the 
regularity of the trajectories based in the vicinity of the corresponding configuration 
(see discussion in section H75)) . This is clearly visible in the leftmost and rightmost 
panels of the figure, pertaining to r = 1.5 and r — 1.82, respectively. 

that shows the Poincare sections obtained from the same set of base points considered 
in the case r = 1.66. Despite some features signaUing a certain degree of chaoticity, the 
sections are much more regular and can be easily distinguished from each other. The 
unstable fixed point, as well as the stable one, disappears also if r is increased by a 
sufficient amount (the disappearance of fixed points takes place at the branch tip after 
a coalescence process). Similar to the previously examined situation, an increase in the 
regularity of the Poincare sections is observed for r = 1.82, when only dimeric fixed 
points are present (rightmost panel of figure [9]). 

5. Phase-triggered population oscillations near the CDW configuration 

The last case considers several unexpected and interesting features of CDW fixed points. 
Remarkably, these fixed points are always present independently from the value of the 
parameters r and z/, and have a simple configuration where the entire boson population 
is equally shared between the lateral condensates, being the phases opposite to each 
other. Figure [5] shows that the CDW are saddle points and are expected to exhibit 
a stable character in extended regions of the parameter plane. This fact, is widely 
confirmed by our numeric simulations where, for initial configurations sufficiently close 
to the CDW fixed point, we observe a regular dynamics exhibiting (possibly quite 
complex) periodic oscillations of the three condensate populations. In particular, the 
lateral condensates cyclically exchange small amounts of bosons, so that their population 
imbalance oscillates. Whereas the central condensate, whose population oscillates 
with a considerably smaller amplitude and a frequency twice than that of the lateral 
condensates, participates in the dynamics acting as a bridge between the lateral ones. 
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Figure HD] shows the evolution of a configuration quite close to the CDW for two 
parameter choices belonging to different stability regions of figure O As expected, in 




150 200 
rescaled time 



Figure 10. Evolution of an initial configuration slightly displaced from the CDW, 
z = (1/2,0, l/2e^^^^'^''), if — 5 ■ 10^'^. Black and grey solid lines represent ni/N and 
n2/N, respectively. Panels a and b refer to the parameter choice v = 0, t — 0.703, 
belonging to a stable region but very close to the border dividing it from an unstable 
region (see figure [5] and (jA.ll[) ). The evolution displayed in panel c was obtained 
increasing r by 0.005. This was sufficient to cross the border and enter the unstable 
region. The predicted instability of the CDW is confirmed by the strongly irregular 
character of the dynamics, which fully shows up after an apparently regular transient. 
The remaining population, n^, (not shown) can be found exploiting the conservation 
of the total number of bosons J^j "-j = the regular case, panels a and b, it 

essentially balances the oscillations of ni, n2 being almost vanishing. 



the stable case (panels a and b), the amplitude of the population oscillations increases 
with increasing displacement of the initial configuration from the CDW. Of course, in 
general, such displacement cannot be increased indefinitely: the initial configuration 
would eventually leave the "domain of infiuence" of the CDW and the consequent 
trajectory would be determined by the stability character of a different fixed point. 
The same considerations apply to configurations obtained from the CDW by setting 
the phase difference between the lateral condensates to a value 631 = \6s — 9i\ 7^ vr and 
maintaining the total boson population equally shared by the lateral condensates. In 
this case the oscillation amplitude is a function of |^3 — with a minimum at |03i| = tt, 
where it vanishes. Our numeric simulations confirm these predictions. In fact, for 
sufficiently small r's, regular oscillatory dynamics is observed over the entire range of 
^31. This is clearly shown in figure fTT| displaying the population oscillations pertaining 
to initial configurations identified by several values of 6*31. 

A further important feature of the CDW fixed points is highlighted by figure [T2l 
showing that for small values of r, the amplitude ai of the oscillation of the population 
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Figure 11. Oscillation of the population ni relevant to initial configurations where 
the total number of bosons is equally shared by the lateral condensates. The three 
panels refer to different (increasing) values of the phase difference 6*31 — {d^ — 9i\. a. 
^31 = 9/107r; b. 831 — l/27r; c. ^31 = I/IOtt. Black and grey hues refer to r = 0.1 
and r = 0.2, respectively. In every case there was no bias on the central condensate 
{v — 0). Notice that the oscillation amplitude increases with increasing 631. Notice 
also the incipient anharmonic shape of the oscillations for small ^31 (panel c.) 



ui exhibits a quite simple dependence on the phase difference ^31 of the relevant initial 
configuration. More precisely ai linearly depends on cos(6'3i/2). A rigorous discussion 
of such dependence is beyond the scope of the present article. However, some insight 
in this sense can be gained by adopting an approach similar to that developed in [17] 
in the study of the stability of the dimeric integrable subregime. Interestingly, the 
same linear behavior as in figure [I2] is obtained by relaxing the requirement that the 
entire boson population is equally shared by the two lateral condensates. This is clearly 
shown in figure [13], where we plotted the average oscillation amplitude relevant to initial 
configurations characterized by a certain phase difference ^31 and condensate populations 
Uj slightly displaced from that of the CDW fixed point. 

The plots displayed in figures [12] and [13] suggest the possibility of measuring the 
phase difference between two equally populated condensates. Let us stressing the fact 
that, in order to profitably exploit this quantum effect to measure the difference of 
phase between two condensates, we need of a nondestructive real-time probing of the 
condensates such as the technique employed to image dimer-like condensates confined 
by magnetic traps [HI [50] . 

6. Conclusions 

The present paper has been inspired by the wish of gaining a deeper insight both 
on the complexity of three-well dynamics (which, as noted above, is the paradigm of 
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Figure 12. Amplitude of the population oscillations of one of the lateral condensates 
for different values of the initial phase difference ^31 and v = 0. Crosses (x): r = 0.1; 
Plus signs (+): r = 0.15. The straight lines are linear fits of the two sets of data. 
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Figure 13. A simple dependence of the population oscillation on the initial phase 
difference 6*31 is present even if the initial population configuration does not reproduce 
that of the CDW exactly. Each of the plotted data was obtained averaging the results 
of 300 simulations where the initial populations were allowed to differ by about 1% 
from the configuration (1/2,0,1/2), and the phase of the central (almost depleted) 
condensate was chosen randomly. 
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the behavior one can expect in more complex BECs arrays) and on the possibihty of 
detecting macroscopic effects amenable to experimental observation. 

Recently, it has been proposed that the superposition of multiple laser beams could 
allow to create unusual optical lattices in a controlled manner ^5T]- In this respect, we 
observe that a suitable choice of laser wavelengths could generate a superlattice whose 
elementary cell is characterized by a three- well potential such as the one addressed in the 
present paper. The further minimization of the tunneling between subsequent trimeric 
cells should provide an alternative experimental realization of the system under concern. 
A nice example in this sense is the trimerized Kagome optical lattice [52j providing 
(multiple copies of) the trapping potential typical of a symmetric closed trimer. 

In the present paper we focus on the rich phenomenology offered by the trimer 
and on the possibility of triggering macroscopic effects through a suitable choice of the 
adjustable parameters, relegating as much as possible to the appendices the technical 
details necessarily involved in a systematic analysis. The macroscopic effects we evidence 
in this paper are obtained through a careful analysis of the stability diagrams introduced 
in [39] , whose operational value is reviewed in section [31 Their derivation is thoroughly 
discussed in Appendix A Such diagrams describe the trimer phase space in terms of the 
stability character of dynamical fixed points for any choice of the significant parameters. 
This wealth of information enriches the recent literature about the dynamical stability of 
mesoscopic BEC arrays [I9]-[2T], [25]- [29], [53] and supplies an essential tool in designing 
experimental realizations of the trimer dynamics. Also, at the theoretical level, it 
provides a sound basis for gaining a deeper insight into the quantum counterpart of 
classical dynamical instability, which is far from being fully understood. 

The use of stability diagrams leads to identify in section H] three different 
mechanisms able to modify in a substantial way the trimer djTiamics. These are 
exploited in sections 14.11 14.21 and 14.31 to prime instability in breather-like regular 
oscillations, loss of coherence in anti-breather states, and suppression of instability due 
to fixed-point disappearance, respectively. Macroscopic dynamical oscillations are shown 
to be triggered by means of slight parameter changes thus providing a stimulating basis 
for planning actual experiments in the scheme described in [l]-[8]. The analysis of the 
stability diagrams in section [5] also evidences that a simple configuration exists allowing 
to relate in a direct manner the initial phase difference between two condensates to the 
amplitude of the ensuing population oscillations. This suggests a rather simple scheme 
aimed at testing experimentally the limits of validity of the mean-field picture, where 
each condensate is endowed with a macroscopic phase. The study of the fixed points 
-inherent in the classical picture- reported in sections [2] and [3] proves useful also in 
unraveling the complex structure of the energy spectrum of the quantum trimer, in that 
it allows to recognize a (either exact or approximate) mapping with a simpler dimeric 
system [5l]. The mapping is shown to persist also at the quantum level, thus extending 
the results found in [2l] to the case of the quantum trimer. These results provide a 
useful starting point in investigating the quantum counterpart [33] of classical dynamical 
instability. A deeper insight into the relation between the classical and quantum picture 
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is certainly useful for the approaches to this issue based on the dynamical algebra method 
[T3| H2I W7\ or on the measure of quantum-state entanglement [55] . 



Appendix A. Derivation of the stability diagrams 



The linear stability character of a fixed point x* = {xi,X2,X3) can be determined by 
studying the eigenvalues of the matrix associated to the (linearized) equations governing 
the dynamics of the (small) displacements v* = (q*, p*) from the fixed point itself, where 
z* = (zi, Z2, Zs) = VN (x + q* + i p*). The Hamiltonian, up to the second order in the 
small displacements v*, reads 



where 



n{z) -xN = n{ 



Mh. 
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and Ai 



2(2x2 



^j2^)/'T- As it is well known, linear instability occurs when 



at least one of such (in principle complex) eigenvalues features a positive real part. 
The linearized equations obtained from Hamiltonian flA.l|) through Poisson's brackets 



2Nh 



(these stem from the original ones {zj^ z^} 



have the form 



UNt 

2h 



M.v, 




(A.3) 



In accordance with the expected time reversal symmetry of the linear dynamics, the 
six (complex) eigenvalues of matrix Ms come in pairs featuring the same modulus and 
opposite signs. More precisely each pair of eigenvalues of Mg consists in the signed square 
roots of an eigenvalue of the three-by-three matrix —QP. Indeed, due to the diagonal 
block structure of Ms, the relevant eigenvalue equation is T>{X) = det(A/ — Ms) = 
det{X^I + QP) = A2(A^ — sA^ + p) = 0, where in writing the characteristic polynomial 
of —QP we took into account the fact that, similar to P, such a matrix has a zero 
eigenvalue. This feature is the consequence of the fact that both q* and p* are 
constrained by the condition -|- |z2p + l^sP = N. According to 'P(A) = the 
fixed point is (linearly) stable provided that the three conditions 



- 4p > 0, 



s < 0, 



p>0 



(A.4) 

are simultaneously met. Indeed in this case the two non-zero eigenvalues of —QP are 
real and non-positive, so that the four non-zero eigenvalues of Mg — namely their signed 
square roots — are purely imaginary. Conversely, whenever an eigenvalue of —QP is 
complex or real positive, one of its square roots has necessarily a positive real part, 
giving rise to (linear) instability. The following subsections describe in some detail the 
construction of the stability diagrams for the DFP, the NFP and the CDW. 
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Appendix A. 1. Dimeric fixed points 

As illustrated in section [2l the dimeric fixed points (DFP) pertain to the integrable 

dimeric subregime [17] of the dynamics described by system ([3]), which is characterized 

by zi{t) = Z3{t). In this situation the fixed point configuration (xi,X2,X3) is described 

by xi(^) = X3{9) = COS0/V2 and X2{9) = sin 9 (see Eqs. (jH])), where the tangent 

a{9) = tan(0) = X2/ {V^ xi) must be a real root of polynomial V{a;T,u) = + 

A^T-^a^ - B^T-^a - 1 (see formula (ED) with A^ = V2{2 - u), = ^/2(1 + z/). 

Before proceeding to a more detailed discussion of the DFP case, we remark that the 

previous polar representation of Xj's allows one to readily relate variable 9 to the fixed 

point {xi, X2, X3). Particularly, tan(6') = X2/{\/2xi) shows how 9^0 corresponds 

to a situation where the central condensate is almost depleted, so that the total 

population is (equally) shared by the lateral condensates. Conversely, when 9 ~ =1=1 =p 

the lateral condensates are almost empty and the total population is concentrated in 

the central one. Notice also that the sign of 9 is similarly strictly related to the phase 

differences between the central condensate and the lateral ones. More precisely such 

phase differences are or ±7r depending on whether 9 is positive or negative. 

As we discussed above, solving the fixed point equations for the DFP amounts to 

finding the real roots of the fourth-degree polynomial V{a; r, u) (henceforth denoted by 

Pd{o:] t]u)). Four real solutions are found from Pd{c(; r; z/) = for 

1 r 1 3/2 

|z/-l/2| < 1/2, T<^ {3/2f^ -\u-l/2f^ (A.5) 

V2 L J 

and two real solutions elsewhere in the parameter plane, as it is illustrated in figure [TJ 
In section El we showed how equation Poid', r; z/) = suggests that, for any choice of 
parameters t = f, u = u, the DFP can be found in the 9 — r plane from the intersection 
of the parametric curve (see [T5!) 

r- (z^ - 2)^2 + (zy + 1) ^ 
— 1 

with the vertical straight line t = f. The curve Tnia; v) features two asymptotes at 
9 = ±7r/4 {a = ±1). This is in agreement with the fact that polynomial PD(a;r, z/) 
has no less than two real solutions. In particular it is easy to verify that when z/ = 1/2 
the values 9 = ±1/2 are solutions of Pd{c(',t]1') = 0. Hence, in this case, the two 
unbounded branches coalesce with their asymptotes, thus providing r-independent 
solutions corresponding to the configurations {xi,X2,xs) = (1/2, ±1/^2, 1/2). More 
generally, when z/ 7^ 1/2 the branches of the curve rD(a;^'), and hence the roots of 
PD(a;; r; u) = 0, are characterized by a reflection symmetry with respect to the nearest 
asymptote (see figure [2]). This can be seen by noticing that 

7r\ 1" 



tan I ?9 ± ^ 



2 z/ T 3 sin(2i9) cos(2i9) ^ ^ 



TT TT 

'4' 4 



(A.7) 



sin(2^?) 2 a/2 

where the auxiliary angle {} describes the positive (negative) solutions 9 = {} + it/4:{9- 
t9 - 7r/4) . One easily verifies that td [tan{^ ± 7r/4) ; z/ + 1/2] = td [tan(±7r /4 - ??) ; 1/2 - z/ 
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Local character of the DFP. This information is obtained by exploiting the 
constraint Xj ^ and the fixed point equations IQ. After expressing the quantities 
appearing in matrices Q and P (see (1A.2I) ) as Ai = A3 = X2/X1 = \pla and 
A2 = [x\ + X'i)/x2 = V^/cy, the analysis of the signature of both matrix Q and P 
depending on parameter a = tan 6 leads to the following conclusions. For any choice of 
the parameters r and z/, a DFP is a minimum if the corresponding angular coordinate, 
6, is positive. If, conversely, < the fixed point is either a maximum or a saddle point 
depending on whether r is greater or less than d{ta.n6) where d{a) = max [di{a), d2{a)] 
with ^1(0;) = — 2V2/[a(l + a^)] and (^2(0) = —6V2/[a^{l + a^)^]. Combining such 
results with the linear stability analysis based on equation T>{X) = (defined after 
formula flA.SP ) confirms that both maxima and minima are, as expected, stable fixed 
points, whereas all the saddles of the dimeric regime are unstable. The situation for the 
DFP is summarized in the corresponding stability diagram displayed in figure [21 The 
differently shaded region pertain to fixed points with a different stability character. 



Appendix A. 2. Non dimeric regime 

Since in this case xi ^ X3, the polar representation of nondimeric fixed points (NFP) 
is not as simple as that of the DFP case. Still, one can take advantage from the 
conservation of the total boson number. Some algebraic manipulations of the fixed- 
point equations (see non- dimeric solutions in section 12. ip suggest to switch to the 
new coordinates Xi = (xi + X3)/-\/2, X2 = X2 and X3 = (xi — xs)/\^, obeying the 
constraint Xf + X| + X| = 1, before setting Xi = RcosO and X2 = RsinO. This way 
X3 = ±a/1 — Xf — xl and radial coordinate (|TT1) are shown to depend on 6 and r as 



, 'l + a2)(2V2-ar) , U + V2a{l + Ta^) , , 

R(a:T) = x- ^ -,X3(a)=±\ — r- (A. 8) 

^ ' ^ ^' 2^2(2 + a2) ' ' y 4(2 + ^ ^ 

with a = tan^ = X2/X1, and the fixed points equations are equivalent to the equation 
for the roots of polynomial V{a; r, u) = a'^ + A^T~^a^ — Bi^r^^a — 4/3 = 0, (henceforth 
denoted by PN{a]u,T)), where A^ = 2^2(2 -z/)/3, = 4^2(1 + z/)/3 (see [I2]). 
Recalling that we are ultimately looking for real solutions of system ([6]) it is clear that 
the roots which make R{a) or X^la) imaginary must be discarded. Before analysing how 
NFP are identified, it is important to remark that in the NFP case the relation among the 
auxiliary variable 6 and the configuration of the fixed point is not as transparent as in the 
DFP case, basically due to the 6 dependence of the radial variable flA.Sp . Nevertheless, 
within the NFP stability diagram some special regions can be still recognized. Notice 
indeed that, according to the previous definitions of Xi, X2 and X3 as well as to (lA.Sp . 



the (forbidden) limit situations X3 and R correspond to the dimeric case 
Xi = Xs and to the CDW case Xi = —x^,, X2 = 0, respectively. As we shall discuss 
below, these limit situations are actually the boundaries of two forbidden regions in 
the NFP diagram (see ( lA.lOp and discussion). Due to the restrictions on R and X^, 



the diagram of figure [3], giving the number of fixed points relevant to system on 
varying parameters r and u, is significantly richer than the corresponding diagram 
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for the number of real roots of polynomial P]\f{a; Ujt). Indeed the latter consists of 
four regions, delimited by the solid lines appearing in figure [31 Two of these regions 
correspond to situations where Pat (a; z/, r) features only two real roots, whereas within 
the remaining regions the roots are either four or none. Despite they are real, some of 
these roots can make either R{a) or X3[a) complex. Hence some new regions appear 
within the ones delimited by the solid curves, where the number of solutions must be 
decreased by either one or two units. 

As discussed in section [3l the features of diagram ([3]) can be better understood by 
resorting to the 9 — t plane, where, similar to the DFP case, for a given choice r = f, 
z/ = z/, the real roots of polynomial Pn{(^'-, ^, t) correspond to the intersections of curve 

r^(a;z/) =2a[(z/-2)a2 + z/ + 2]/(3a^+ 1), z/ = z> , (A.9) 

with the straight line r = f. In the half plane r > such curve consists in general 
of two bounded, roughly bell-shaped, branches. More precisely they are defined within 
the 9 intervals [— z/fe,0] and [z/b,7r/2], respectively, where Vh = arctan(A/ (2 + z/)/(2 — z/)). 
When z/ > 2 (z/ < —2) the former (latter) branch, which from now on will be referred 
to as the negative (positive) branch, disappears. These features readily explain the 
four regions in the diagram describing the number of real roots of Pn{o:] J^, t). Hence, 
provided f is sufficiently small, there are either four or two solutions depending on 
whether |z>| < 2 or |z>| > 2. Indeed, for any straight line at a given value of r, there 
exist for each branch a value of z/ such that the branch and the straight line are tangent. 
Denoting such values i^^ir), where the superscript labels the branch, it is easy to check 
that ^ufi^r) is always a growing function of r. Hence two solutions relevant to the 
negative (positive) branch are always present for z/ > z/j~(t) [u < z/^^(r)]. Notice further 
that the presence of a region of parameters where four solutions are found implies that 
there exists some Tc > such that i^ti'^c) = ^ti'^c)- Then it is clear that PN{a; Ujt) 
features four real roots in the region {UfT < u < v^] r < Tc}, no real roots in the 
region {vf < u < uf; r > rd and two real roots elsewhere. The explicit expressions of 
functions z/j^(t) and z/j^(t) for a < and a > 0, respectively, (solid lines in figure HI) can 
be easily described in a parametric way by means of two curves Tj(q;), z/j(a): the curve 
Tt{a) = 16a^/(3a;^ + 9a^ — Sa"^ — 1) intersects TN^a] y) at the tips of its bell-shaped 
branches while z/j(a) is such that ^TAr(a; z/) = at z/ = Vtio) and rAr[a; z/i(a;)] = r((a). 
The requirement that i?, X3 G M determines two forbidden regions in the stability 
diagram:{6' < 0, r > r3(tan0)} and {0 > 0, r > r^(tan6')}, where 

r3(a) = -2/a(l + 2a^), rjf(a) = 2/a, (A.IO) 

(see figure H]). The information about the forbidden regions can be visualized in parallel 
in the t — v plane by drawing the curves z/j:j(r) and ^^{t) describing the intersections 
between rjv(a; v) and the borders of the forbidden regions, tr^q) and T^{a), respectively 
(dotted curves of figure [3]). According to what we said so far, these curves enclose 
two regions where one of the real roots of polynomial PN{oi]i'^T) must be excluded 
from the set of solutions of system ([6]): v > v^ij) and z/ < z//j(r). Similarly two of 
these real roots must be discarded in the regions {Tij(z/) < r < r^^y), v < v^} and 
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{ts^u) < t < (z/), z/ > z/j }. Notice also that tj:j(z/) [t3(z/)] is expectedly tangent to 
^H^) Ki^)] at r = T+, [t = T+]. 

Local character of the NFP. The study of the signature of both matrix Q and P 
defined in ([QD reveals that NFPs are either maxima or saddle points depending on 
whether they lie within the region {9 < 0, ri(tan6') < r < r3(tan6')} or elsewhere 
in the allowed zone. In parallel, the linear stability analysis based on polynomial 
T>{\) (see the discussion relevant to flA.4p ) discloses the (linear) stability of most 
of the saddle points lying in the positive 6 quadrant, namely the ones belonging to 
the light grey region of figure HI Indeed, within the allowed region {6 < 0, t < 
r3(tan^)} U {9 > 0, t < Tji{ta,n6)}, the stability condition (1A.4I) fails to apply in 
the region {6 < 0; r < Tt{ta,n6)} U {6 > 0; t > Tt{tan9)}, where p < 0, and in a 
peculiarly shaped region of the positive 6 quadrant, where s'^ — 4p < 0. Since s is 
positive everywhere outside these regions, it turns out that linear stability characterizes 
not only, as expected, the maxima, but also most of the saddle points at positive 9. 



Appendix A. 3. Central depleted well 

As discussed in section [21 there is only one parameter- independent CDW fixed point: 
X = (1/v^, 0, — 1/v^), which, according to the signature of matrices Q and P (see 
IA.2I) . is always a saddle point. Indeed, no matter for the sign of A2 = —2 (1 + z/)/r, 
the quadratic form has always an undefined signature. Recalling that Ai = A3 = 0, 
the coefficients involved stability condition (1A.4I) have the quite simple form s = 4[r^ + 
(1 + z/)2]/r2 and p = A[t^ + 4(1 + iy)]/T^ yielding A = 16 [2 (z/^ - 1^ + {u + 1)^] /t\ 
Hence, introducing the functions 

ci(z/) = 2 V^T^, C2(z/) = v/(l + ^^)V[2 (1 - z^)] (A.ll) 

the CDW saddle is actually stable in the parameter regions {r > ci(z/); z/ < — 1}, 
{r < C2(z/); |z/| < 1} and {t > 0; u > 1}, and unstable elsewhere. Notice that the 
former zone can be also described as z/ < z//j(r), where z//j is strictly related to the 
boundary of one of the forbidden zones of the NFP stability diagram. Actually, as we 
already noticed above, the points lying on such (forbidden) boundary r = Tji{tan6) 
have a CDW form, since R{a; r) = there (see (1A.10|) and the defining formulas for 



Xj = Xj{R,6), j = 1,2,3 in section Appendix A.2). The stability diagram for the 



CDW fixed point is displayed in figure [51 Notice that, unlike figures [2l and [H it is in 
the plane of the parameters, r and z/. This is possible because, as we already remarked, 
there is always a single, parameter-independent CDW configuration. 
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